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We have studied the nucleation and the physical proper- 
ties of a -1/2 wedge disclination line near the free surface of a 
confined nematic liquid crystal. The position of the disclina- 
tion line has been related to the material parameters (elastic 
constants, anchoring energy and favored anchoring angle of 
the molecules at the free surface). The use of a planar model 
for the structure of the director field (whose predictions have 
been contrasted to those of a fully three-dimensional model) 
has allowed us to relate the experimentally observed position 
of the disclination line to the relevant properties of the liquid 
crystals. In particular, we have been able to observe the col- 
lapse of the disclination line due to a temperature-induced an- 
choring angle transition, which has allowed us to rule out the 
presence of a real disclination line near the nematic/isotropic 
front in directional growth experiments. 
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I. INTRODUCTION 

Disclination lines jjj are characteristic defects of uni- 
axial nematic phases. These lines break the orientational 
order of the director field. They can easily form, for in- 
stance by stirring a nematic drop deposited on a glass 
plate. Observation in the microscope reveals a threaded 
texture composed of thin and thick threads which are 
primarily ±1/2 and ±1 wedge disclination lines. The 
former have singular cores whereas the latter form con- 
tinuous configurations since the director can escape along 
the disclination axis (the so-called escape in the third di- 
mension [0). These lines cost an energy per unit length 
of the order of K (the relevant Frank constant), so they 
tend to spontaneously collapse when the flow vanishes. 
They can even all disappear at the end if the constraints 
at the boundaries and in the bulk are negligible (no ex- 
ternal field). By contrast, some of them may subsist if 
the boundary conditions impose topological constraints 
which are incompatible with a continuous director field. 
Such a phenomenom occurs when the liquid crystal is 
confined inside a sphere with strong homeotropic an- 
choring (nematic molecules anchored parallel to the local 
normal of the surface). In this case, a point defect oc- 
curs at the center of the sphere (hedgehog configuration) 
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HQ . This configuration is stable if the anchoring energy 
per unit surface area, W, (see below for a more detailed 
definition of this quantity) is strong enough to prevent 
the homeotropic anchoring to be broken. More precisely, 
a defect will appear if the anchoring penetration length 
K/W [[| is smaller than the sphere radius R. Such con- 
figurations can be obtained by dispersing a nematic phase 
in a non miscible isotropic liquid || . They are also found 
in polymer dispersed liquid crystals (PDLC) 0. These 
materials, of great practical interest, are obtained by a 
polymerization-induced phase separation with an aver- 
age droplet radius typically ranging between 0.1 fan and 
10 /im. Another way to nucleate a disclination line is to 
confine the nematic liquid crystal inside a capillary tube 
with a strong homeotropic anchoring. If the tube is cir- 
cular, the line is a +1 disclination [||. On the other hand, 
two +1/2 disclinations must occur near both sides of a 
flat rectangular capillary tube ||. In all these examples 
the nematic phase forms a simply connected medium. 
Conversely, it is possible to disperse a non-miscible liq- 
uid into a nematic liquid crystal and form this way an 
emulsion. This problem has recently attracted attention 
JlOt because the elastic distorsions produced by the point 
defects that nucleate near each droplet create long range 
dipolar forces between the droplets. 

In this article, we have chosen to study the nucleation 
and the physical properties of the -1/2 wedge disclination 
line that forms at the edge of a nematic sample confined 
between two parallel glass plates treated chemically to in- 
duce a homeotropic anchoring. At the edge, the nematic 
sample is in contact with an isotropic fluid. The interface 
between the nematic phase and the surrounding isotropic 
fluid is generally curved across the sample thickness be- 
cause the two phases wet the glass plates differently. The 
first experimentally addressed example was that of the 
nematic interface in contact with its own melt. Such a 
situation is achieved in directional solidification when a 
homeotopic nematic sample is placed in a temperature 
gradient so that the nematic/isotropic phase transition 
temperature lies somewhere in the sample. Previous ex- 
periments in directional solidification have shown that a 
disclination line can detach from the interface in the form 
of triangular domains |ll]] . To explain this observation, 
it was suggested that a disclination line lies behind the 
meniscus separating the two phases and that a dust par- 
ticle crossing the line could detach it from the meniscus, 
leading to the observed triangles. This interpretation was 
reinforced by the fact that energy considerations alone 
suggest this line can exist near the meniscus. Another 
reason to suppose its existence was that the wavelength 
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measured at the onset of the Mullins-Sekerka instability 
was five times larger than the wavelength predicted by 
the linear stability analysis ]l2|-|l4[ , suggesting a possible 
renormalization of the apparent surface tension due to 
the disclination line and to the associated elastic effects 
|l5| , [l6f . Unfortunately the line was not visible in the mi- 
croscope in the thin samples used in directional growth. 
For this reason we tried to see the line by increasing the 
sample thickness (its distance to the meniscus should es- 
sentially scale with the sample thickness). Surprisingly, 
even though we increased the sample thickness up to 400 
/xm, we never saw the line and concluded that it might 
be virtual. 

Prompted by this failure, we attempted to observe a 
disclination line at a free nematic interface (i.e., in con- 
tact with air), and we discovered that the line may or may 
not exist in that system, depending on both the strength 
of the anchoring energy and, much more important, on 
the value of the preferred tilt angle of the molecules at 
the free interface. 

This article is structured as follows: in section (||) 
we detail the basic experimental setup and procedure, 
together with the main phenomenology involved in the 
formation of a disclination line near a free interface in 



a confined nematic phase. In section (HI) we outline a 
planar model for the director field in the presence of a 
disclinat ion lin e which we use to do a numerical study. In 
section ( III C ) a more realistic, three-dimensional struc- 
ture of the director field is considered and we show it is a 
good approximation to use the planar model to calculate 
the position of the disclination line. Finally, in section 
( |lV|) we consider the effect of the anchoring angle on the 
disclination line, and use our planar model to study a 
temperature-induced anchoring angle transition. 



II. EXPERIMENT 

We have prepared thin nematic (N) liquid crystal sam- 
ples confined between two parallel glass plates. In order 
to favor a homeotropic anchoring of the liquid crystal 
on the glass plates, thus inducing a homeotropic align- 
ment of the liquid crystal in the bulk, we have treated 
the surface of the glass plates by chemically adsorbing a 
surface active layer on them |17| . Prior to the chemical 
treatment, the glass plate is thoroughly cleaned by gently 
washing it with soap and water in order to remove any 
dirt spots, followed by soaking in soapy water for about 
60 minutes, then rinsed in distilled water and quickly 
dried, first with clean air, and then baked under vacuum 
at 110°C for about 30 minutes. We have tested the use 
of CTAB (hexadecyl trimcthyl ammonium bromide) and 
the use of a silane compound from Merck (ZLI-2510) to 
align the liquid crystal molecules. In order to treat the 
plates with CTAB, we dip them into a 0.002% aqueous 
solution of the surface- active agent. The plates are subse- 
quently lifted slowly out of the solution in order to avoid 




FIG. 1. Cross section of the translationally invariant dis- 
tortion caused on the nematic director field by a boundary 
with an isotropic medium. The distortion may generate a 
— 7r-disclination line along the meniscus (D). The presence of 
the isotropic medium induces the liquid crystal molecules to 
be anchored on the meniscus at a given angle ((f) with respect 
to the local normal (y). 8 is the angle between the normal and 
the horizontal axis. We assume that the meniscus has a circu- 
lar profile with radius R = 6/(2cos(a)), and a wetting angle 
on the glass walls a. The director field distribution is given, 
in the planar approximation, by n = (cos($), 0, sin(<I>)). The 
origin of coordinates is taken to be at the tip of the N/air 



meniscus. 



precipitation of CTAB. For the treatment with the silane 
compound, we prepare a saturated solution of the surface 
active agent in toluene. The glass plates are dipped in the 
solution for a few seconds and quickly lifted out to avoid 
creating dry silane spots. The excess silane is removed 
from the glass plates by gently rubbing with lens clean- 
ing tissue. The cells are mounted with the desired gap 
between plates using either nylon thread (for gaps in the 
range 100 /im - 400 /im) or nickel wire (for gaps smaller 
than 100 /im) as spacer. A precision better that 5% in the 
cell gap is obtained this way. The cells are partially filled 
by capillarity, keeping them at a temperature in the ne- 
matic range of the liquid crystal. Since the observations 
are to be made in the nematic phase, we make sure that 
the sample does not go into the isotropic phase once the 
cells are filled. Even though homeotropic alignment of 
the nematic phase is achieved faster if the product in the 
cell is first melted into the isotropic phase, we have ob- 
served in some cases that the motion of the N / air menis- 
cus caused by the melting of the nematic phase crucially 
alters the observable phenomena. Sealing of the cell can 
also be performed to avoid degradation of the material 
or to prevent motion of the N/air meniscus. Different 
liquid crystal mesophases have been used in our experi- 
ments: 5CB (K15 from Merck), 8CB (K24 from Merck), 
MBBA (from ACROS), and 152 (from Merck). Samples 
with gaps in the range 10/xm < b < 450/im have been 
used. 
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FIG. 2. Disclination line observed from above near a 8CB 
N/air interface confined between two parallel glass plates. 
The cell gap is b — 100 fim. The sample is illuminated with 
polarized light. H: homeotropically aligned nematic phase. 
P: planar aligned nematic phase between the disclination line 
and the meniscus. M: N/air meniscus from whose width we 
can estimate the wetting angle on the glass plates (see Fig. |^). 
We can observe a disclination line that is parallel to the menis- 
cus, except for the region near a three-dimensional structure 
(pointed to by an arrow). The presence of such defects sug- 
gests that the director field has a structure that escapes the 
plane formed by the normal to the meniscus and the perpen- 
dicular to the glass plates. 



A. Preliminary observations 

When a disclination line is present, it appears equidis- 
tant from the two glass plate boundaries and parallel to 
the meniscus, suggesting a structure that is translation- 
ally invariant along the interface (see Fig. [l]). When the 
sample is observed between crossed polarizers (with the 
polarizers being parallel or perpendicular to the discli- 
nation line), the region between the disclination line and 
the homeotropic bulk reveals the presence of a distor- 
tion of the director field which indicates that it has a 
component in the direction of the disclination line (see 
Fig. H). The presence of three-dimensional defects on 
the disclination line, and the fact that the position of 
the line is altered in the vicinity of such defects suggest 
that the three-dimensional structure of the director field 
may play a significant role in the position of the discli- 
nation line. As we will see below, however, a simple 
planar model for the disclination line (in which the di- 
rector field has no component along the line) is able to 
reproduce both qualitatively and quantitatively our ex- 
perimental observations on the line position (see section 
HI). A three-dimensional extension of the planar model 
shows that even though the component of the director 
field along the disclination line is not negligible, its ef- 
fects on the posi tion of the line are not very significant 
(see section III C ). 

In order to measure the distance of the disclination line 
to the meniscus, we always select a region that is reason- 
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FIG. 3. Distance of the disclination line to the meniscus in 
a 5CB N/air interface. The distance is proportional to the 
gap between the glass plates, and d/b — 0.094 ± 0.005. The 
inset shows how the wetting contact angle of the meniscus 
at the glass walls is experimentally measured by assuming 
a circular profile and by measuring the width of the menis- 
cus, w (see Fig. |l|). An average value of w/b is obtained 
using measurements with all cell gaps. This average value 
(slope of the solid line) is used to compute the contact angle, 
a = 7r/2 — 2tan _1 (2w/6) (a ~ 25° for this configuration). 



ably far from any such defects and that does not appear 
distorted by them. For gaps b < 10 /im, both the limit 
of our optical resolution and the presence of the three- 
dimensional defects make our observations difficult. We 
are able to observe the formation of the disclination line, 
but we are unable to meaningfully measure its distance 
to the meniscus. Therefore, we will restrict our measure- 
ments to samples with b > 20 /im. 



B. Equilibrium position of the disclination line 

We first studied the dependence of the position of 
the disclination line on the thickness of the sample. 
The measurements are performed at a constant tem- 
perature in a region of the nematic phase where tem- 
perature variations of about 1°C have a negligible ef- 
fect on the position of the disclination line. Typically, 
this is the situation when the temperature is more that 
1°C above the Nematic/Smectic-A phase transition (pro- 
vided a Smectic-A phase exists) and more that 1°C be- 
low the Nematic/Isotropic phase transition. In Fig. || 
we plot the sample-thickness dependence of the position 
of the disclination line for 5CB. Our liquid crystal has 
TV/ — 34.7°C, and we perform our measurements at 
T ~ 30.0°C (T NI - T ~ 4.7°C). As we see in Fig. |, the 
distance of the disclination line to the meniscus scales 
linearly with the cell gap over the range of gaps we have 
explored. Clearly, the shape of the meniscus, particu- 
larly the wetting angle with the glass plates, will have a 
strong influence on the position of the disclination line. 
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TABLE I. Ratio of distance of disclination line to N/air 
meniscus, d, to the cell gap, b, and wetting contact angle 
on the glass plates, a, for several liquid crystals and surface 
treatments. The wetting properties of the glass plates are 
noticeably different with the two different surface treatments 
we have used. 



ZLI-2510 


5CB 8CB 


MBBA 152 


d/b 
a 


0.094±0.005 0.105±0.005 
25° 20° 


0.f05±0.005 0.096±0.005 

25° 8° 




CTAB 


8CB 


MBBA 


d/b 
a 


0.100±0.005 
f0° 


0.100±0.005 
f0° 



Assuming a constant radius of curvature (circular pro- 
file), we are able to experimentally estimate the wetting 
angle on the plates (see Fig. H). Given that the N/air 
surface energy (~ 10~ 3 — 10~ Jm~ 2 ) is orders of mag- 
nitude larger than the nematic anchoring energy on the 
free interface (W ~ 10~ 6 — 10~ 5 Jm~ 2 ) |jL8| , a shape with 
a circular profile, which both minimizes the surface area 
between the two media and meets the glass plates at a 
given contact angle, should be expected. Similar results 
are observed with 8CB, MBBA, and 152, with a combina- 
tion of d/b ratios and wetting contact angles a as shown 
in Table ffl. It is customary to construct the dimension- 
less parameter Wb/K, which compares the length scale 
imposed by the confined geometry (in our case the cell 
gap, &), with an intrinsic length scale of the material, 
K/W. Note that throughout this article we will suppose 
that the anchoring energy on the glass plates (which in- 
duces the homeotropic orientation of the bulk) is large 
compared to W, so that the presence of the free interface 
does not alter the anchoring condition on the plates. 

The perturbation caused on the homeotropic direc- 
tor field configuration by the presence of the free in- 
terface will be less apparent for samples thinner than 
K/W. One expects the formation of a disclination line 
for Wb/K » I, when surface anchoring effects domi- 
nate and impose a rigid anchoring angle. If Wb/K << 1, 
however, elastic effects will dominate and a balance be- 
tween elastic and surface anchoring torques, rather that 
a rigid anchoring angle, is expected to determine the 
anchoring condition at the meniscus. In this case, the 
distortion caused by the free interface will be less im- 
portant, and the equilibrium structure may be one with- 
out a disclination line. The results reported in Fig. || 
and Table |, namely, a disclination line whose position 
is completely determined by the geometry of the system, 
prove that the only relevant length scale is b, the cell 
gap. We are clearly in a regime with Wb/K >> 1, or 
infinite anchoring energy, where the anchoring angle on 
the meniscus is rigidly set. If we were able to make clear 
observations with a small enough sample thickness (typ- 
ically b < 10 /im), we should presumably see the end of 
the d (x b regime, and observe the collapse of the discli- 
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FIG. 4. Collapse of the disclination line close to Tna for 
8CB in the presence of a N/air interface with a sample of 
thickness b = 25 /im. The imposed temperature gradient is 
38.6°Ccm _1 . Top: picture of the sample from above, ob- 
served between partially crossed polarizers. The length of the 
vertical line segment is 10 and that of the horizontal seg- 
ment is 15 fim. Tna occurs at the left edge of the picture. 
Bottom: the distance of the disclination line to the meniscus 
is measured as a function of the distance to the transition 
point, and converted to T — Tna by means of the known tem- 
perature gradient. The disclination line reaches its plateau 
position at about 0.5 K above Tna- 



nation line. 

In the next section, we outline a model to relate the 
equilibrium position of the disclination line to the dif- 
ferent experimental parameters. Our model predicts 
the disclination line to start collapsing noticeably for 
Wb/K < 5 (see Fig. g), which we may use to give a 
lower bound for the anchoring energy, W. We have con- 
fidently measured the position of the disclination line 
for b > 25 /im, and still observe the same d/b ra- 
tio (within error bars). We can, thus, estimate W > 
hK/b = 1.1 10" 6 Jnr 2 {K 3 = 5.5 10" 12 N for 5CB at 
Tni — T — 5.0°C |l9|]). A similar estimate is obtained 
from our measurements with 8CB. 

In an attempt to observe experimentally the collapse 
of the disclination line due to small values of Wb/K, 
we have considered altering the value of W by replacing 
the air interface with an interface with another isotropic 
medium, such as an isotropic fluid or glass wires with 
different surface treatments. Even though we have ob- 
served an effect of the different interfaces on the anchor- 
ing conditions, the configurations had poor reproducibil- 
ity and dubious stability. We also tried to increase Wb/K 
by approaching the nematic/Smectic-A phase transition 
temperature. Indeed, we know from previous work that 
the bend elastic constant K$ diverges at this transition 
p9| fH[ so that the disclination line energy also diverges. 
The experiment was performed by placing a 25 /im-thick 
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sample of 8CB inside a constant temperature gradient 
parallel to the liquid-crystal/ air interface (see Fig. ||). 
The immediate observation is that while in the N phase 
a disclination line is observed, it disappears as we en- 
ter the SmA phase. More precisely, d/b shifts from its 
asymptotic value (close to 0.1) when T — T^a < 0.3°C 
and typically equals 0.05 at T - T NA ~ 0.1°C. At this 
temperature K3 ~ 4 10~ n N while d/b = 0.05 when 
Wb/K ~ 1 (see the following section). This allows us to 
roughly estimate W in 8CB: W ~ 1.5 10~ 6 Jm~ 2 . 



III. MODEL 

In order to better relate the phenomena described 
in the previous section to the experimental parameters, 
we modeled the director-field configuration in a con- 
fined nematic geometry in the presence of a free in- 
terface. Even though our observations show a three- 
dimensional director-field structure around the disclina- 



tion line (see section II A), we propose a simple, two- 
dimensional model with a planar disclination line (the 
director field is assumed to be perpendicular to the line). 
This model has allowed us to obtain numerical estimates 
for the position of the disclination line that match very 
well the experimental observations. Finally, we con- 
sider a more accurate model, where the three-dimensional 
structure of the director field is taken into account. We 
verify that, in fact, the predictions from the planar model 
quite accurately predict the position of the disclination 
line. 
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dx I \ dz 



(2) 



where K\ ~ K 3 = K. The glass walls are taken to be at 
z = ±b/2, and we assume that a strong anchoring fixes 
$ = ±7r/2 there. 

The interaction with the N/air meniscus contributes a 
surface energy that is customarily approximated by 



W 2, x 

7 = 7o - y cos (ip- ifo), 



(3) 



with ip = <f> — 8. Here, 9 is the angle between the lo- 
cal normal to the meniscus and the horizontal direction, 
and ipa introduces a favored anchoring angle between the 
director field and the local normal to the interface for 
which the surface energy density is minimum (see Fig. 

The energy introduced in the system by the presence 
of the N/air interface will be 

E =\ K j (grad$) 2 dv-^W J cos 2 (<^ - po) dS, (4) 

where the volume integral spans over the bulk nematic 
phase and the surface integral spans over the N/air 
meniscus. 

We will seek the director-field distribution (character- 
ized by $) by minimizing the energy functional with re- 
spect to Integrating by parts and imposing that the 
director field be homeotropic both on the glass bound- 
aries and far from the interface, one obtains 



A. Equations 

If we consider that the distortion caused by the pres- 
ence of the free interface is translationally invariant along 
the interface, and ignore the component of the direc- 
tor field along the disclination line (planar disclination 
line), then a typical director field configuration will ap- 
pear as shown in Fig. |], and it will be characterized 
by the angular field $(x, z), with n = (cos($), 0, sin(<I>)). 
Given the two-dimensional nature of this approximation, 
there will be no contribution from the twist of the di- 
rector field to the elastic free energy (n • curl n = 0). 
Likewise, the saddle-splay term (K24) will vanish iden- 
tically 0,0. Moreover, for the experimental conditions 
described in the previous section in which the samples 
are kept at a temperature away from both the N/I and 
N/SmA transitions, the splay and bend elastic constants 
take on roughly the same values, Ki ~ K3. A final ap- 
proximation is made in which the mixed splay-bend term 
(Ki 3 ) is neglected 22 2§]. With all this, we can write, 



A$ = 0, in the bulk, 



(5a) 



/ = -K [(div n) 2 + (n x curl n) 2 ] , (1) 



<9$ W 

K-qZ — — sin[2(</? — V?o)]- on the free interface. (5b) 



Eq. (5b) establishes a torque balance at the inter- 
face (here, d/dv represents the derivative along the local 
normal direction, and 6 is the angle between the local 
normal to the interface (u) and the horizontal direction). 
The shape of the interface is completely determined by 
the experimentally measured wetting contact angle (see 
Fig. ||). These equations can be made dimensionless by 
rescaling all length scales by b and defining the anchoring 
parameter Wb/K, which contains all the material prop- 
erties. With this, one would expect the distance from 
the disclination line to the meniscus to have a functional 
dependence of the form d = bf(Wb/K). In the strong an- 
choring regime (Wb/K — > 00), <f> = 8 + cpo will be rigidly 
set on the interface, and d oc b is expected. Note that 
there strictly exists in our model a discontinuity in the 
director field at the tip of the meniscus when ipo > is 
rigidly imposed at the interface. This may lead to values 
of d smaller than the real ones. The use of eq. 5b with a 



which gives 



large, but finite, value of Wb/K effectively removes this 
discontinuity and shows that the error introduced by the 
infinite anchoring approximation is very small. 
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FIG. 5. Expected position of the planar disclination line, 
d, in terms of the different experimental parameters and ex- 
pressed as the ratio d/b. (•): fixing a = 25° and ipo = 0°, the 
position of the disclination line as a function of the anchoring 
strength parameter, Wb/K, is calculated, d/b reaches a lim- 
iting value as Wb/K is increased (limit of infinite anchoring). 
In the insets we show the ratio d/b (in the Wb/K — > oo limit) 
as a function of the wetting angle, a (□), in the ipo = con- 
figuration, and as a function of the anchoring angle, ipo (A), 
for a wetting angle a — 25°. Note that our model predicts 
that the disclination line will approach the meniscus roughly 
linearly with increasing ipo, and that one can extrapolate the 
value of ipo at which the line would collapse (dashed line). 



B. Numerical model 

In order to find the configuration which minimizes the 
total energy, we proceed to solve eq. || by relaxing the 
fixed boundary problem obtained by setting the position 
of the disclination line at a distance d from the menis- 
cus. Therefore, 4>(x, O*) = ±7r/2 is set for x > d, while 
$(x,0) = is set for < x < d. In order to avoid the 
7r-discontinuity at z = 0, we restrict our calculation to 
the z > region. The value of the wetting angle (which 
can be determined experimentally, see Fig. g| and Table |) 
and the value of the favored anchoring angle, ifo, together 
with the condition of homeotropic configuration both far 
from the meniscus and on the glass walls, completely de- 
termine the director field distribution. Typically, we use 
discretization grids of about 400x200 points. Once the 
configuration is obtained, we compute its total energy 
and then minimize this energy with respect to the config- 
urations obtained by changing the value of d. Note that 
we neglect the contribution of the core of the disclination 
to the total energy. We assume that it will be constant, 
independent of d and, therefore, will not be relevant in 
the determination of the configuration with the least en- 
ergy. In Fig. ||, we show the results of these calculations. 
First, we explore d/b as a function of Wb/K. This corre- 
sponds to the measurements performed in section ( II B ). 
We observe that the model predicts that d/b converges 



asymptotically to d/b = 0.105 ± 0.003 as Wb/K -> oo, 
for a = 25°, and tpo = 0. Note that this is, within er- 
ror bars, consistent with the results reported in Table 
|| Only our measurements with 5CB (shown in Fig. |^) 
are slightly smaller that the prediction from our model. 
Our numerical calculations predict as well that, as Wb/K 
gets smaller, d/b starts to decrease. In order for d/b to 
decrease more than, say, 10% below the infinite anchor- 
ing position (which would be experimentally measurable 
with our optical resolution) one would need Wb/K < 5 
(see Fig. ||) . We have also explored the effect that a and 
ipo have on the disclination line. We can see in Fig. ^ 
that d/b increases with increasing values of a. Indeed, 
if we held d/b fixed, a flatter meniscus (larger a) would 
require larger distortions of the director field in order to 
match both the strong anchoring on the glass walls and 
on the meniscus. This would increase the elastic energy, 
which can be prevented by increasing d. The experimen- 
tal observations reported in Table | are consistent with 
this result, even though the dispersion in the data, com- 
bined with the fact that in the a < 25° region d/b is 
not very sensitive to changes in a, do not allow for any 
conclusive trend to be extracted from our measurements. 
On the other hand, we may exploit the modeled effects of 
the anchoring angle, ipo, to study the known phenomenon 
of temperature-induced anchoring transitions (see section 
IV). We observe (see Fig. ||) that, within our numerical 
resolution, d decreases roughly linearly with increasing 
values of (fo- Using this result, one can immediately cor- 
relate changes in the position of the disclination line with 
changes in the anchoring condition, provided W remains 
unchanged. Observe that, extrapolating the calculations 
reported in Fig. |5|, we can expect that, even in the pres- 
ence of a strong anchoring, a configuration without a 
disclination line may be observed, provided ipo is large 
enough (ipo > 25° for the configuration in Fig. |J). This 
is the reason why a real disclination line is not formed 
spontaneously near a nematic/isotropic front, since on 
that interface large values for ipo are reported (<po — 48° 
for 8CB Hj). We recall that the anchoring length at this 
interface K/W is too small (of the order of 5fim with 
W ~ 8.5 10~ 7 Jin" 2 H and K 3 ~ 4 10~ 12 N (l9))to ex- 
plain alone the absence of a real disclination line in thick 
samples (d ~ 50 — 400/im). We also emphasize that the 
formation of planar triangles in directional solidification 
PJ requires that a dust particle crosses the N-I interface, 
suggesting the disclination line is not real but virtual (its 
core lies in the isotropic liquid). This observation fully 
agrees with our predictions. 



C. Three dimensional effects 



As we have seen in section (II A), our experimental 
observations reveal a three-dimensional structure in the 
director field in the vicinity of the disclination line. The 
model described above assumes a planar configuration 
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IV. EFFECT OF THE ANCHORING ANGLE 




We have seen on section [I] that, based on the results 



FIG. 6. Three dimensional structure of the director field 
near the N/air meniscus. The field lines are tangent to the 
planar component of the director field, and the shading rep- 
resents the angle between the director field and the plane, O, 
where = is the planar case. Darker areas correspond to 
larger third-dimensional component. The same angle is plot- 
ted as a surface plot inside the box. The escaped component 
diverges in the core of the disclination line. The core of the 
disclination is centered at a distance d ~ 0.116, in agreement 
with the planar model and with our observations. 



of the director field, in which the component along the 
disclination line is ignored. In this section, we study the 
modifications introduced by a three-dimensional model 
of the director field. 

We have used the full Landau-Ginzburg-de Gennes ten- 
sorial expansion for the elastic free energy density in 
power series of the quadrupolar order parameter (which 
we define as Qij = (Q/2)(3njnj — <%), assuming a uni- 
axial system). We take K\ = K% = 2K2, which is an ex- 
cellent approximation for 5CB and 8CB in the region of 
the nematic phase we study. Moreover, the saddle-splay 
term is included as well (K24, = —K\) and the mixed 
splay-bend term is set to zero — 0) 23 1. We have 

only considered the case of rigid homeotropic anchoring, 
in which the liquid crystal molecules are forced to anchor 
on the meniscus at cpo = 0° . A detailed description of the 
model and the numerical methods involved are given in 
©■ 

Indeed, we obtain an equilibrium configuration where 
the director field has a significant projection in the di- 
rection of the disclination line, especially in the vicinity 
of the line (see Fig. ^) . Surprisingly, though, when the 
position of the disclination line is estimated, the result 
is, within error bars, the same as the one found with 
the much simpler planar model. Therefore, the three- 
dimensional calculations validate the use of the planar 
model to predict the position of the disclination line as a 
function of the relevant physical parameters. 



from our computations, we expect the disclination line 
to approach the meniscus as the anchoring angle is larger 
than zero and that, presumably, configurations without 
a disclination line might be possible, provided the an- 
choring angle is large enough. It is well-known that the 
equilibrium anchoring angle of the nematic molecules at 
a N/air interface can depend on temperature [p6|-p8| . We 
have performed measurements with MBBA and 152, both 
of which are known to have such a behavior, in order to 
study the effect the anchoring angle has on the formation 
of a disclination line. Samples are prepared as explained 
above and placed in a temperature-regulated oven and 
the position of the disclination line as a function of tem- 
perature is measured. 



A. Continuous anchoring angle transition in MBBA 
and 152 

For 152, the disclination line starts to approach conti- 
nously the N/air meniscus as we increase the tempera- 
ture beyond To — 25° C, until it reaches a final configu- 
ration where the distance remains unchanged by further 
increasing the temperature. No hysteresis is observed as 
we alternate heating and cooling cycles, and the tran- 
sition temperature is independent of sample thickness. 
In fact, measurements performed with different values of 
the cell gap are collapsed by rescaling the distance of the 
disclination line with the cell gap (see Fig. (7J) . This sug- 
gests that what we are seeing is not due to changes in the 
anchoring energy, since its influence in the position of the 
disclination line is introduced by the parameter Wb/K, 
which depends on the sample thickness. For MBBA, the 
distance from the disclination line to the N/air meniscus 
starts to decrease continuously as we lower the tempera- 
ture below To ~ 26° C. Just as with 152, no hysteresis is 
observed here. The transition temperature is slightly dif- 
ferent for different samples, however, and it may depend 
on sample thickness, even though we have been unable to 
find any meaningful trend. Shifting the temperature axis 
with the value of To for each experiment, and rescaling 
the position of the disclination line by b effectively brings 
onto a single curve the data taken with different sample 
thicknesses (see Fig. 0) suggesting as well that we are 
always in the infinite anchoring regime. 

Both for 152 and for MBBA, we see that the depar- 
ture from the upper plateau is well-described by a con- 
tinuous (or supercritical) bifurcation of the form d/b oc 
|T— Tol 1 / 2 , which is analogous to what has been observed 
elsewhere for anchoring-angle transitions p8| . Fitting 
such an analytical profile to our data, we can get a more 
precise estimate for the transition temperature, T). 

Note that, even though a second stable position for 
the disclination line is observed for 152, the line keeps ap- 
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FIG. 7. Distance of the disclination line to the N/air 
meniscus (d), rescaled by the sample thickness (b) as a func- 
tion of the departure from To, the temperature above or be- 
low which the disclination line starts approaching the menis- 
cus, showing a continuous transition between two different 
configurations presumably due to an anchoring angle tran- 
sition. Top: measurements performed with 152, with sam- 
ple thicknesses b — 50 /im, 110 /jm, 140 /jm, and 220 /im. 
Bottom: measurements performed with MBBA, with sample 
thicknesses b = 110 /im, 190 /im, and 220 /im. The departure 
from the ipo = configuration (d/b ~ 0.1) is compatible with 
d/bcx \T - To | 1/2 (solid curve). The solid horizontal lines are 
guides for the eye in the plateau regions, where d/b does not 
change with T. In both systems, heating and cooling cycles 
are indistinguishable, as long as the disclination line does not 
collapse (MBBA case). The diagram in the low temperature 
region illustrates the hysteretical collapse of the disclination 
line observed for this product. 



^ a ^ 



FIG. 8. A nematic microdroplet of 152 is placed on a sub- 
strate that induces a homeotropic anchoring and it is ob- 
served between crossed polarizers, at different temperatures: 
T = 34.5°C (a), T = 28.2°C (b), T = 26.2°C (c), and 
T = 25.2° C (d). The temperature of (a) is well in the lower 
plateau on Fig. 0. The observation of a point defect on the 
free surface suggests a tilted anchoring of the molecules. As 
the temperature is reduced, the distortion decreases until we 
lower T below To ~ 25.5° C so that we are in the upper plateau 
on Fig. [?]. The distortion has vanished, suggesting we have a 
homeotropic anchoring at the free surface. 



proaching the meniscus as the MBBA samples are cooled 
to lower temperatures, until the disclination line finally 
collapses. As we describe below, there is hysteresis in this 
collapse, and the dynamics that leads to it is reminiscent 
of a first-order transition. Based on the observations re- 
ported above and on the results found in the literature, 
we would like to argue that the changes we see in the 
position of the disclination line when we vary the tem- 
perature of the samples can be accounted for by changes 
in the anchoring angle at the interface. The existence of 
an anchoring angle transition is easily demonstrated by 
observing between crossed polarizers a micro-droplet of 
the liquid crystal placed on a substrate that induces a 
homeotropic anchoring (see Fig. |J). For temperatures 
at which the disclination line is at the upper plateau in 
Fig. 0, we have observed no point defects in the droplets. 
As we vary the temperature entering the region where a 
transition in the position of the disclination is observed 
in Fig. 0, a characteristic point defect structure starts to 
become manifest, and the distortion in the director field 
increases in the temperature region where the disclina- 
tion line is closest to the meniscus. We have used our 
planar model for the structure of the director field near 
the N/air interface in order to relate the position of the 
disclination line to the anchoring angle (see Fig. |). The 
numerical model considers a planar disclination line and 
a strong anchoring condition at the meniscus with an an- 
choring angle that can be different from the homeotropic 
orientation (ipo = 0). Note that our numerical model 
assumes K\ ~ K$. For MBBA, in the relevant temper- 
ature range (T ~ 25° C), one finds in the literature that 
0.84 < K3/K1 < 1.4. The same information, however, 
is not available for 152. The fact that our planar model 
successfully predicts the position of the disclination line 
in the homeotropic anchoring situation is a good indica- 
tion that the above approximation is valid. Moreover, in 
the vicinity of the continuous anchoring angle transition, 
the angular dependence of the surface energy may have 
its leading order term cancelled. Indeed, the surface en- 
ergy in the vicinity of the second order anchoring angle 
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FIG. 9. Temperature-induced orientational transition of 
the anchoring angle at the N/air interface where the angle 
has been estimated by transforming the d/b vs. T — To data 
on Fig. with our numerical model. Top: results corre- 
sponding to the data obtained with 152 (To ~ 25° C) and 
measurements reported by Sonin et al. [H (solid line). Bot- 
tom: results corresponding to the data obtained with MBBA 
(□) (To ~ 26° C), results reported by Sonin et al. (•), and 
results reported by Bouchiat et al. |27j (solid line). 

transition can be expanded as 

7M=7o + -^(T-T V 2 + iV + ..., (6) 

where Wq and B are two positive constants. The 
equilibrium angle is <po = for T > Tq and po = 
^W {T -T )/B for T < T (MBBA case). The torque 
originated by an anchoring at an angle tp > at T = Tq 
is 

p- = W (T-T )<p + 4Btp 3 + .... (7) 

We see that, at T = To, the leading order term vanishes. 
We will assume that the first non-vanishing term, 4Btp 3 
still increases steeply when ip > 0, so that the anchoring 
angle remains close to its value at equilibrium. 

The numerical calculations show that the dependence 
of the ratio d/b on the anchoring angle can be roughly 
approximated by a straight line, whose slope depends on 
the wetting angle of the meniscus on the glass plates (see 
Fig. ||). Using our numerical results, we have been able 
to transform the d/b vs. T — To data in Fig. ^| into 
ipo vs. T — T (see Fig. ||). Among other features, our 



po vs. T — To data maintains a continuous bifurcation 
compatible with tp oc |T — TqI 1 / 2 , consistent with pre- 
vious measurements by other groups. Our estimations 
for the dependence of the anchoring angle on tempera- 
ture for MBBA are well in agreement with all the results 
reported in the literature, both at the onset of the bi- 
furcation and at much lower temperatures (^7), ex- 
tending almost 10°C below the transition. We observe a 
homeotropic anchoring above To and a steadily increas- 
ing anchoring angle whose saturation we are unable to 
observe as we enter a regime where the disclination line 
collapses. We see that, provided the anchoring angle is 
large enough ((po > 25°) a configuration without a discli- 
nation line is possible. On the other hand, our estimates 
for 152 depart from the homeotropic anchoring configu- 
ration for T above To, with values for ipo that are about 
1.7 times lower than what was reported by Sonin et al. 
pil. For T well above To, we find a stable configuration 
with ipo ~ 16° tilted from the homeotropic orientation. 
The limited available information on the physical prop- 
erties of 152 may account, at least partially, for these dis- 
crepancies, even though our only assumptions are strong 
anchoring and one-constant approximation for the elastic 
energy. 

B. Collapse of the disclination line in MBBA: a 
first-order transition 

When the temperature is lowered below a thickness- 
dependent value (typically 10 degrees below T ), regions 
where the disclination line has collapsed start to nucle- 
ate around dust particles on the meniscus (see Fig. fio| ). 
These collapsed regions propagate until they eventually 
extend along all the interface, and the equilibrium con- 
figuration has no disclination line. Before collapsing, the 
disclination line in the now metastable regions keeps ap- 
proaching the meniscus as the temperature is lowered 
even further. We have already explained that, as long as 
the disclination line does not collapse, there is no hys- 
teresis between heating and cooling cycles. Once the line 
has collapsed, however, heating the sample does not re- 
generate the line continuously. Instead, there is a hetero- 
geneous nucleation process similar to the one that leads 
to the collapse of the line upon cooling, and the temper- 
ature at which the line is regenerated is well above the 
temperature at which the line started to collapse. The 
diagram in Fig. (?] illustrates this hysteresis cycle. 

This dynamics suggests the coexistence of two possi- 
ble configurations in the structure of the director field 
for MBBA, one presenting a disclination line and another 
with a continuous distortion of the director field. The ob- 
served nucleation mechanism for the collapse and regen- 
eration of the disclination line shows that the transition 
between these two configuration is energy activated. A 
detailed study of this process would yield a good quanti- 
tative characterization of the energy associated with the 
core of the disclination line, which has yet to be done. 
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FIG. 10. Hysteresis observed in the nucleation and in the 
collapse of a disclination line near the N/air interface for a 
6 = 220 /im MBBA sample. The cycle is illustrated in the di- 
agram in Fig. (?]. Top (al-a3) : collapse of the disclination line 
after a temperature quench of a few tenths of a degree into a 
final temperature of T ~ 16.3° C. The configuration contain- 
ing a disclination line becomes unstable, and regions without 
a disclination line nucleate and propagate. The collapsing 
disclination line is pinned on a dust particle. Time sequence: 
t = (al), t = 8 s (a2), t = 16 s (a3). Bottom (bl-b5): nu- 
cleation of the disclination line observed upon heating of the 
same sample, once the disclination line had totally collapsed. 
The disclination line starts to regenerate around a dust parti- 
cle at T ~ 19.3° C and it propagates while T is held constant, 
filling the whole interface. The pinning defect on the disclina- 
tion line marks where the dust particle is. The time sequence 
is as follows: t = (bl), t = 48s (b2), t = 96s (b3), t = 144s 
(b4), t = 192 s (b5). 

V. SUMMARY AND CONCLUSIONS 

In this article, we have studied the nucleation and the 
physical properties of a -1/2 wedge disclination line which 
forms near the open edge of a nematic liquid crystal 
sample confined between two glass plates that induce a 
homeotropic orientation of the nematic phase. The po- 
sition of the disclination line, in particular its distance 
to the edge of the N/air meniscus is directly measur- 
able, and we have related it to the physical properties 
of the liquid crystal with the help of a planar model 
for the structure of the director field in the presence 
of the N/air interface. Since the three-dimensional na- 
ture of the structure of the director field in the vicinity 
of the disclination line is experimentally evidenced, we 
have used a three-dimensional model to obtain a more 
accurate description of this structure. Indeed, we have 
observed a significant escaped component of the direc- 



tor field, even though the position of the disclination line 
that minimizes the total energy agrees with the predic- 
tions of the planar model, which we have therefore used 
since it allows to implement realistic boundary conditions 
in a simple way. We have performed different attempts at 
altering the position of the disclination line by changing 
the anchoring parameter Wb/K. For large values of this 
parameter, the position of the disclination line is com- 
pletely determined by the geometry of the problem, and 
therefore d oc b. By changing the sample thickness, we 
have tried to reduce the magnitude of Wb/K so that the 
line aproaches the meniscus faster than the linear depen- 
dence on b. Our attempts at observing this regime and 
extracting a value for W this way have failed, most likely 
due to the experimental limitation in the smallest values 
for b we have been able to use. Under usual experimental 
conditions, the disclination line will form and its position 
will agree with the predictions from the infinite anchor- 
ing approximation. Altering W in a way that can affect 
the position of the disclination line has proven to be a 
difficult task. Given the fact that W is an intrinsic prop- 
erty of the interface, one would need to consider different 
isotropic media or even chemically altering the anchoring. 
We have been unable to alter the nature of the anchor- 
ing strength in a stable and reproducible way. Another 
promising way to control the magnitude of Wb/K is ap- 
proaching the Nematic/Smectic-A phase transition tem- 
perature. The bend elastic constant, K3, can be made 
arbitrarily large by approaching that temperature. By 
placing our samples in a temperature gradient, we have 
been able to observe the collapse of the disclination line, 
which we see to start approaching the N/air meniscus 
when the temperature is less than a few hundred mil- 
likelvin above Tjva- From this experiment, we estimated 
W in 8CB near Tna- Nevertheless, this method has im- 
portant experimental limitations (presence of a vertical 
temperature gradient in the samples, existence of a zig- 
zag instability of the disclination line at very small tem- 
perature gradient, bad knowledge of the elastic constants 
near the transition, etc.). Finally, our calculations show 
that a departure of the anchoring condition at the N/air 
meniscus from the homeotropic orientation will signifi- 
cantly alter the position of the disclination line. Indeed, 
we have been able to observe the disclination line to ap- 
proach the meniscus when materials known to have a 
temperature-induced anchoring angle transition are used. 
Moreover, we have shown that a configuration without 
a disclination line is obtained for values of the anchor- 
ing angle above 25° (a value that depends slightly on 
the material properties), even though energy considera- 
tions alone would suggest the possibility of the forma- 
tion of a disclination line (Wb/K is large). We believe 
this is the situation we encounter in the vicinity of a ne- 
matic/isotropic liquid crystal front in directional growth 
experiments. Even though we are in a regime of infinite 
anchoring (large Wb/K), a large value of the anchoring 
angle (about 45°) guarantees that the disclination line 
will be virtual in those cases. 
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